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To date, several experimental groups reported measurements of the thickness dependence of Tc 
of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the Tc-oscillations 
varies significantly from one experiment to another. Here we propose that the reason for this 
unresolved issue is an interplay of the quantum-size variations in the single-electron density of 
states with thickness-dependent oscillations in the phonon mediated electron-electron coupling. Such 
oscillations in the coupling depend on the substrate material, the quality of the interface, the 
protection cover and other details of the fabrication process, changing from one experiment to 
another. This explains why the available data do not exhibit one-voice consistency about the 
amplitude of the Tc-oscillations. Our analyses are based on a numerical solution of the Bogoliubov- 
de Gennes equations for a superconducting slab. 

PACS numbers: 74.78.Na 



I. INTRODUCTION 

Quantum-size effects on the excitation gapi^ and the 
transition temperature^"— in superconducting ultrathin 
films have been a subject of study since Blatt and 
Thompson in 1963 predicted that the superconducting 
gap in ultrathin films oscillates as a function of the thick- 
nessji Such oscillations are caused by the quantization 
of the electron motion in the direction normal to the 
film. In particular, the formation of discrete electron 
levels, i.e., the quantum well states (QWS))ii^ splits the 
three-dimensional (3D) conduction band into a series of 
two-dimensional (2D) subbands. The bottom of each 
subband is at the energy position of the corresponding 
QWS and so it is shifted down (up) when increasing (de- 
creasing) the nanofilm thickness. When one of the QWS 
approaches the Fermi level, the density of the single- 
electron states at the Fermi level iV(0) increases, resulting 
in enhanced superconductivity. 

Recently, superconducting Pb nanofilms have been 
intensively investigated experimentally where electronic 
effects are much stronger than effects due to strain. 
The critical temperature of atomically flat Pb nanofllms 
with thickness down to a few monolayers on a Si(lll) 
substrate have been measured under different condi- 
tions)^"— i^"— It was found that superconductivity is not 
destroyed by fluctuations even when the thickness of the 
nanofilm is only a single monolayer^ and the depen- 
dence of the critical temperature on thickness shows os- 
cillations with a period close to two monolayers (ML),— 
i.e., the even-odd staggering in the superconducting prop- 
erties. 

In addition to N{0), the phonon-mediated coupling g 
between electrons also makes an important contribution 
to superconductivity. Large effortaii"— have been made 
to understand how the formation of QWS can affect the 
electron mass-enhancement factor A ^ gN{0) [only in 
the weak-coupling limit A = giV(O)]. Experimentally, 
from the information of the temperature dependence of 



the line widths of QWS in Ag films deposited on top of 
Fe(lll) whisker, it was found that A exhibits an over- 
all decrease with increasing film thickness so that A ap- 
proaches its bulk value as « l/N, where N is the number 
of monolayersiii The experimental study of A in Pb(lll) 
nanofilms on silicon showed that the mass-enhancement 
factor increases with increasing TV. In both cases A ex- 
hibits pronounced thickness-dependent oscillations. As it 
is mentioned above, a well-known reason for these oscilla- 
tions is the variation of N{Q) with the nanofilm thickness. 
However, this is not the only important reason. There is 
an additional (less known) effect coming from thickness- 
dependent oscillations of the coupling g. The interplay 
of the oscillations of g and N{0) is not well understood 
yet and often overlooked. In our previous paper— we ini- 
tiated a study of such an interplay based on intuitive ar- 
guments and used an ansatz for the thickness dependent 
g. In the present paper we revisit the problem and de- 
velop a simple model to justify this ansatz and to analyze 
all the currently available experimental data of Tc versus 
film thickness. Our investigation is based on a numerical 
solution of the Bogoliubov-dc Gennes (BdG) equations. 
However, instead of solving the full BdG formalism,— we 
make use of the Anderson recipe for a very accurate ap- 
proximate semi-analytical solution to the BdG equations: 
corrections to the Anderson approximation are less than 
a few percent for nanofilms. This will significantly sim- 
plify our numerical analysis. 

Our paper is organized as follows. In Sec. [11] we out- 
line the BdG equations for a superconducting slab, in- 
cluding the basic moments concerning the Anderson ap- 
proximation. In Sec, mil we discuss our numerical results 
on the critical temperature in ultrathin superconducting 
nanofilms and highlight important aspects due to the in- 
terplay of the thickness-dependent oscillations of A^(0) 
and g. We demonstrate that this interplay is respon- 
sible for significant variations in the amplitude of the 
quantum-size oscillations of Tc from one experiment to 
another. Our conclusions arc given in Sec. IIVI 
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II. FORMALISM 

A. Bogolibov-de Gennes equations for a 
superconducting slab 

The translational invariance is broken in the pres- 
ence of quantum confinement so that the superconduct- 
ing order parameter A becomes position dependent in 
nanoscale superconductors, i.e., A = A(r). In this case 
the profile of A(r) can be found from the BdG equations 
(see the textbooki^) which are written as (for zero mag- 
netic field) 



Air) 
A(r) -i?,, 



■"ik(r) 



(1) 



where = — 2m — /i is the single-electron Hamilto- 
nian (for zero magnetic field) measured from the chemi- 
cal potential ^; i is the quantum number associated with 
QWS for electron motion in the z direction perpendicular 
to the film, and k = (kx, ky) is the wavevector for quasi- 
free motion of electrons parallel to the film {k = |k|); 
Uik(r) and Vik(r) are the particle- and hole- like wave 
functions; Eik stands for the quasi-particlc energy; and 
A(r) is the superconducting order parameter (chosen as 
real). The BdG equations are solved in a self-consistent 
manner with the self-consistency condition given bj*i^ 

A(r) = 5^u,k(r)<k(r)tanh(/3i?,,/2), (2) 

ik 

where g > is the effective coupling of electrons and (3 = 
l/iksT). The summation in Eq. ^ runs over states with 
the single-electron energy ^^fe inside the Debye window, 
i.e., I^ifcl < hwD, with hjOD = 8.27 meV for Pb (see the 
textbooks -i^^) 

As known, the BdG equations can, in principle, be 
solved without additional approximations. However, 
even for ultrathin nanofilms this can be a time-consuming 
task, especially in the vicinity of Tc, where a large number 
of iterations is usually needed for a proper convergence. 
This is why it is useful to invoke Anderson's recipe for a 
semi-analytical approximate solution. The use of Ander- 
son's recipe is well justified because (i) corrections were 
found to be less than a few percent for nanofilms (see, 
e.g., Ref. [ish and (ii) resulting equations arc much easier 
for numerical implementation. In terms of the particle- 
and hole-like amplitudes, the Anderson approximation 
reads 



Uik(r) = ZijfeVik(r), Wik(r) = VifcV'ik(r), 



(3) 



where lAik and Vik are multiplicative factors (real) 
and '!/'ik(r) is the single-electron wave function, i.e., 
-ffeV'ik(r) = CifcV'ik(r). In the case of interest, i.e., for 
a superconducting slab, we have 



^,^^(r) ^ ^^^^^(k^x+kyy) ^ 



(4) 



where (pi (z) is the QWS wave function and Lx , Ly are 
the dimensions of the unit cell used for the periodic 
boundary conditions in the x and y directions. For infi- 
nite confinement of electrons in the slab with thickness 
d {Lx,Ly 3> d) we have ipi{z) = yj2/ ds\\\{-K(i -\- l)z/d) 
and 
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Inserting Eq. ^ into Eq. ([T]) and introducing the 
subband-dependent gap [in the case of interest A(r) = 
A(z)] 



A, = dz^Uz)Mz), 



(6) 



we find a homogeneous system of two linear equations 
controlling Uik and Vik- The determinant of the corre- 
sponding matrix should be equal to zero for a nontriv- 
ial solution, leading to Eik = \/ S^fk + ■ Then, by us- 
ing the standard normalization condition Uff. + Vfj^ = 1, 
the BdG equations are reduced to the BCS-like self- 
consistency equation 



with the interaction-matrix element 



= .g dzipf,{z)iff{z). 



(7) 



(8) 



As T approaches the critical temperature Tc, Eik re- 
duces to the single-electron energy £^ik independent of 
Ai. Hence, at Tc Eq. ([7]) becomes a linear homogeneous 
system of equations for A.^ that has a nontrivial solution 
when the corresponding determinant is zero, i.e.. 



det[M,,i-Si,i] 
from which we obtain Tc- In Eq. 



0, 



(9) 



I de0(e-£o||tanh(/3,|^|/2), (10) 



with /3c — l/{kBTc) and Si 



^(^+1) 



fi. We note 
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that a numerical solution of Eq. (jH]) is significantly faster 
than the procedure based on a numerical study of the 
BdG equations. The point is that Eq. ([9]) yields directly 
the value of Tc- In contrast, to calculate Tc from the 
BdG equations, one first needs to find A^ as function 
of temperature. Then, Tc can be extracted from these 
data as the temperature above which there exists only 
zero solution for A^. In addition, such a BdG-based nu- 
merical procedure of determining Tc is rather capricious 
in practice because at temperatures close to Tc conver- 
gence becomes rather slow and it is difficult to estimate 
an appropriate number of iterations to approach a solu- 
tion. This often results in an overestimation of Tc that is 
larger than the corrections to the Anderson approxima- 
tion, i.e., one more solid argument in favor of the present 
numerical scheme. 
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B. Quantum-size variations of the electron-electron 
coupling 



In Sec. Ill Al it is assumed that the phonon mediated 
couphng of electrons g is not position dependent. It is 
true for bulk and a good approximation even for metal- 
lic nanofilms with thickness d > 10-20 nm (see our re- 
sults below). However, this is not the case for ultrathin 
nanofilms due to significantly different (as compared to 
bulk) conditions for the lattice vibrations at the inter- 
face between the film and the semiconductor substrate. 
As we are interested in superconducting quantum-size 
oscillations typical for ultrathin metallic films, the ques- 
tion arises how our formalism should be modified to take 
account of this feature. 

The coupling between electrons will depend on the 
proximity to the interface, i.e., g = g{z). Deep in the 
film, g{z) approaches the bulk coupling g^. For Pb we 
takei^ 5oA^buik(0) = 0.39, with A^buik(O) the bulk den- 
sity of states at the Fermi energy. Approaching the in- 
terface, g{z) is not equal to g^ any longer but acquires a 
different value ga. Therefore, to first approximation we 
can replace g{z) by the step function 



.gif < z < dif, 
go dif < z< d. 



(11) 



where djf can be interpreted as the interface thickness. 
When introducing the spatial dependence of the coupling 
between electrons, Eq. ([5|) is replaced by 



so that we obtain 



<izgiz)(pl{z)ip'^{z), 



(1 + '^.'.), 



dzg{z). 



(12) 



(13) 



As seen, to take into account the interface effect on the 
electron-electron effective coupling, one should replace g 
in all the relevant formulas in Sec. Ill Al by its spatially- 
averaged value g given by Eq. ([H]). Based on Eqs. pT|) 
and (dH), we can find 



g = .90 



c^if(.9o - ffif) 
aN 



(14) 



where a is the lattice constant in the direction of the 
film growth and N is the number of monolayers in the 
film {d = aN, and wc take a 0.286A for Pb(lll), 
see Ref. [T2h . As seen from Eq. (fT4)) . the presence of the 
interface results in an extra contribution to g decaying 
with increasing thickness as 1/iV. This is typical of the 
surface effects whose relative contribution can be esti- 
mated as roughly proportional to the ratio of the whole 
surface area to the system volume. Such a contribution 
to the coupling was first reported in Ref. [ll|, where the 
electron-phonon mass enhancement parameter A in Ag 



films grown on top of Fc(lOO) was experimentally studied 
and found to dependent on the film thickness. As known, 
A = gN{0) for g -> 0, with iV(0) the density of states, 
and A is often called the (dimensionless) coupling con- 
stant. Ab initio DFT-calculations of A also revealed the 
presence of the 1/A^-contribution, see, e.g., Refs. fill and 
[Til . Note that the well-known softening of the phonon 
modes due to surface effects (see, e.g., Ref. [l^) can be 
incorporated in a similar manner. Below we assume that 
all surface effects are included in g^f. The same holds for 
effects due to a protective cover. 

While go is the film-material constant, ga is not and 
changes from one experiment to another. Moreover, gif 
varies also with changing nanofilm thickness because of 
quantum-size oscillations in the interface properties, see, 
for instance, theoretical results and discussion about the 
surface energy and the work function in nanofilms.— 
Generally, the quantum-size oscillations are controlled 
by a multiple-subband structure appearing in the single- 
electron spectrum due to the formation of QWS. Each 
time when the bottom of a subband (QWS) crosses fi, a 
shape (quantum-size) resonance appears^ with a signifi- 
cant effect on the system characteristics. From Eq. ([S]) 
it is clear that the shape resonances are periodic with 
period half the Fermi wavelength A_f/2. Thus, assuming 
that (7if exhibits thickness-dependent oscillations with pe- 
riod Ai?/2, Eq. (|14p can be represented in the form 
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gijAiraN/XF) 
N 



(15) 



where gi{x) is a periodic function with gi{0) = (7i(27r). 
Below, in our numerical study of the BdG equations, 
we replace g in the formulas of Sec. Ill Al by g given by 
Eq. (113. 

It is worth noting that one should distinguish oscilla- 
tions of g from those of the electron-phonon mass en- 
hancement parameter A. As A sa gN{0), there are two 
sources of quantum-size oscillations in A, i.e., (i) the den- 
sity of states that exhibits an increase each time when a 
new QWS crosses fi, and (ii) the phonon-mediated cou- 
pling g. The oscillations of A are a known point that was 
investigated previously in experiments and in theoretical 
studies, see, e.g., Refs. [lll-[l4l However, effects of the in- 
terplay of oscillations of A^(0) with those of g were not 
investigated and reported. 

Experimentally measured values of the period of the 
quantum-size oscillations in Pb nanofilms are found to be 
very close to 2 ML. However, the period is not exactly 
2 ML so that the even-odd staggering in the basic su- 
perconducting properties has phase-slip points where the 
phase of the even-odd oscillations changes by tt. When 
4a/Ai? = 1, i.e., the period is equal to 2 ML, then we 
deal only with two parameters 5i(7r) and (7i(27r) be- 
cause gi{ATraN/Xp) is equal to .gi(7r) for an odd num- 
ber of monolayers and to 171 (27r) for an even number of 
ML's. On the contrary, when Aa/Xp is not an integer, 
we need to go into more detail about the dependence 
of gi on N. There is the possibility to overcome this 
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FIG. 1. (Color online) Experimental data on the thickness 
dependence of the critical temperature Tc of Pb nanofilms on 
Si(lll) from different experiments. 



problem by using, say, a "minimal scenario" when choos- 
ing the chemical potential such that the period of the 
quantum-size oscillations is exactly 2 ML. In this case 
we are able to limit ourselves to only two fitting pa- 
rameters .91 (tt) and 5i(27r). Therefore, our theoretical 
analysis will be restricted to the data between the phase- 
slip points where the effects of the mismatch between 
Af/2 and 2 ML can be neglected. Within such a mini- 
mal scenario, the parabolic band approximation used in 
Eq. ([5]) dictates that fi = 1.150 eV. It corresponds to A_f 
being four times the single monolayer thickness a, i.e., 
Xf = 1.14 nm. 

Now, we have everything at our disposal to check if 
the interplay between thickness-dependent oscillations of 
A'^(O) and g can provide a solid understanding of the sig- 
nificant variations in the amplitude of the quantum-size 
oscillations of Tc found in different experiments. 



III. RESULTS AND DISCUSSIONS 

A. Brief review of experimental results on the 
Tc-oscillations 

In Fig. [T] we show all reported sets of recent exper- 
imental data on the thickness dependence of Tc in ul- 
trathin Pb(lll) nanofilms on silicon. The data indi- 
cated by circles are from Rcf. and were obtained by ex 
situ resistivity measurements. Here Pb nanofilms were 
deposited on a Si(lll) substrate (with the Si(lll) (7 x 
7)— reconstruction) and covered by a protective Au layer 
with thickness 4 ML. This set of data exhibits clear 
even-odd quantum-size staggering for thicknesses of 22- 
28 ML. For smaller thicknesses there are only data for 
odd number of atomic layers, i.e., TV = 15, 17, 19 and 



21. A clear phase slip of tt in the even-odd oscilla- 
tions can be observed at = 22 ML. At this point the 
even-odd staggering in Tc changes its trend: above the 
thickness N = 22 ML the critical temperature increases 
when passing from an odd number of monolayers to an 
even number, however it decreases from A = 21 ML to 
N = 22 ML. Below we will limit the analysis of this 
data-set to A > 23. 

The results of Ref. are given by stars and obtained 
from ex situ magnetic-susceptibility measurements for Pb 
films on Si(lll) with the interface Si(lll)(A/3 x V3). 
The samples were protected by a Ge cap. The set in- 
cludes data for films with an odd number of atomic layers 
A = 5, 7, 9, 11, 13 and with even numbers A = 14, 16 
and 18. Here it is not possible to obtain significant in- 
formation about the even-odd oscillations of Tc'. we can 
only compare the data for A = 13 and A = 14. As seen 
from this comparison, the increase of Tc from A = 13 
to A = 14 is about two times smaller than the averaged 
amplitude of the even-odd oscillations observed in the 
previous set of data from Rcf. iSj. 

The third experimental set of thickness dependent Tc 
(rhombuses, Rcf. @) is for Pb nanofilms on the substrate 
Si(lll)(7 X 7). Here the data were measured by in situ 
scanning tunneling spectroscopy. There was no protec- 
tive layer covering the Pb film in this case. Clear sig- 
natures of even-odd quantum-size oscillations are visible 
down to thickness of 5 ML. However, the amplitude of 
these oscillations is about 2-3 times smaller than that of 
the first data-set)^ It is worth noting that Tc of the third 
data set does not show an overall decrease with decreas- 
ing film thickness. It is seen that there are two points 
where the phase of the even-odd staggering of Tc slips by 
TT, i.e., at A = 6 and at A = 14. For this set of data we 
will focus on the domain 7 < A < 13. 

The fourth (squares) and fifth (triangles) data sets are 
from the same group*^ and were measured by in situ 
scanning tunneling spectroscopy. They correspond to 
Pb nanofilms on a Si(lll) substrate with reconstruc- 
tion Si(lll)(7 X 7) and Si(lll)(y3 x ^3), respectively. 
The fourth set shows the even-odd staggering of Tc for 
A = 14-18. However, here the slip-point is situated at 
A = 16, which significantly reduces the number of points 
available for our analysis (below we work with data for 
A < 15). The fifth experimental set for Tc docs not 
provide us with a clear pattern of even-odd oscillations. 
Only from the data for film thicknesses A = 5, 6 and 7 
we can guess that the amplitude of such oscillations is 
not pronounced. The phase-slip point of the even-odd 
staggering is visible at A = 6 (below we make a numer- 
ical fit to the data for A > 7). We note that the data 
of the fourth set are in significant variance with those of 
the third set in spite of the fact that the type of mea- 
surement and the substrate-film interface were similar in 
both experiments. 
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FIG. 2. (Color online) Theoretical results for the even-odd oscillations of Tc (open symbols) versus the experimental data (solid 
symbols): (a) theoretical data were calculated for giiji) = l.Qlgo (odd numbers) and g\(2-K) = 2.13t;o, the experimental set is 
from Ref.H (b) theoretical data were calculated for 31 (tt) = 0.64gio (odd) and gi{2ix) = lAQgo (even), experimental results are 
from Ref. |j; (c) the same as in the previous panels but for (?i(7r) = — 0.12(;eff and gi{2n) = 0.84(;cff (where (;e//Afbuik(0) = 0.36) 
and the experimental set is from Ref. 0; (d) the same as in previous panels but for gi(7r) = 1.08(/o and (/i(27r) = 1.54(/o and the 
experimental set is from Ref.[^ (for the interface Si(lll)(7x 7) — Pb); (e) the theoretical data were calculated for gi(7r) = l.OSpo 
and (7i(27r) = 1.59(?o and the experimental set is from Ref. (for the interface Si(lll)(\/3 x y/i) — Pb). 



B. Impact of the thickness-dependent coupling 
between electrons 



Results of our numerical analysis (open symbols) are 
given versus the experimental data (solid symbols) in 
Fig. [5] In panel (a) we compare our data with the Tc- 
valucs from the first set (Ref. [^. Here we choose the 
range N = 23-28, i.e., above N = 22, the experimental 
phase-slip point of the even-odd oscillations in Tc (see 
the discussion in the previous section). Our calculations 
have been performed for 51 (tt) = 1.67.go (odd numbers) 
and gi(27r) = 2.13go (even numbers), which yields good 
agreement with experiment. At first sight these values of 
(71 seem too large to produce a deviation of only about 
10%-20% from bulk. However, gi/N is the correction 
to the electron-electron coupling rather than gi . In par- 
ticular, for iV = 20 we obtain gi/N = 0.084,9o for odd 
numbers and gi/N = O.llgo for even numbers. It is 
worth noting that gi is larger for even numbers and so 
the resulting coupling go — gi/N is smaller. However, 
Tc is generally larger for even numbers as compared to 
that for odd numbers in Fig. [Ija). The reason for this 
counterintuitive behavior is that there is an essential in- 
terplay between the even-odd oscillations in the density 
of states N{0) and similar oscillations in the coupling 
g. In the case of panel (a) the drop in g for even N is 
fully overcome by an increase in iV(0) so that the prod- 
uct gN{0) [Tc - e-i/(s^(°»] slightly increases for even- 
number films and decreases for odd-number films. 



Figure [21b) shows our theoretical results versus the ex- 
perimental data of the second set from Ref. 0. Here we 
used the fitting parameters 51 (tt) = O.6450 (for odd num- 
bers) and .gi(27r) = 1.46go (for even numbers). As seen, 
gi for odd numbers is three times smaller than for even 
numbers but Tc does not exhibit any pronounced even- 
odd staggering. The explanation is again the interplay 
of the quantum-size variations of iV(0) with thickness- 
dependent oscillations in g: the effect of the difference 
between gi{n) and gi{2'K) is almost compensated due to 
significant drops of N{0) for even numbers of atomic lay- 
ers. It is worth noting the surprisingly good agreement 
over the full range of experimental data. 

Now we switch to the analysis of the third set of experi- 
mental data for Tc from Ref. y. These data are compared 
with our theoretical results in Fig. ^c) for iV = 7 ML 
to TV = 13 ML (i.e., between the two phase-slip points 

= 6 and A^ = 14). For this data-set Tc does not ex- 
hibit an overall decrease with decreasing N but oscillates 
around a value that is about 15% lower than the bulk crit- 
ical temperature of 7.2 K. Even for A^ — 50 ML the criti- 
cal temperature was found to be about 15% smaller than 
in bulk. The bulk value of Tc and the superconducting 
energy gap were achieved only for relatively thick films 
with A^ = 500 ML. To incorporate this feature, we intro- 
duce a slightly smaller effective electron-phonon coupling 
constant g^ff which replaces the bulk coupling constant 
go- By choosing 5offA^buik(0) = 0.36, 5ii(7r) = -0.125off 
and (7i(27r) = O-SAg^s, we compare our results (the open 
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symbols) with the experimental critical temperature data 
(the solid symbols) . As follows from Fig. ^c) , we again 
obtain very good agreement with the experimental data. 
As .91 (tt) < 0, the total electron-phonon coupling con- 
stant gcs — gi/N for the odd-layered nanofilms is slightly 
above the effective-bulk value 0.36 and increases when de- 
creasing N. This is totally different from our analysis of 
the previous data sets, and the direct consequence of this 
feature is that Tc does not show an overall decrease with 
decreasing nanofilm thickness. As seen, in this case the 
thickness-dependent oscillations of the electron-electron 
coupling "kills" such a decrease typical of the remaining 
experimental data sets. Note that the appearance of the 
positive interface-induced contribution to g for the odd- 
layered films, is not a surprise. For example, a similar 
effect was found for Ag films deposited on top of Fc(lOO) 
whisker. 

Figures W[d) and (e) show the comparison of our nu- 
merical data with data (the fourth and fifth sets, respec- 
tively) for Tc from Ref. |9|. For Pb nanofilms grown epi- 
taxially on Si(lll)(7 x 7) [panel (d)] we performed cal- 
culations with the fitting parameters giij^) = l.OSga and 
gi{2n) = 1.54^0- For films with interface Si(lll)(\/3 x 
•\/3) — Pb we used giiir) — l.OSgo and 5i(27r) = 1.59go- 
For the fourth set we focus on the range below the phase 
slip point TV = 16 [N < 15). For the fifth set we work 
above the phase-slip thickness = 6 {N > 7). Here Tc 
increases for the even-layered nanofilms and decreases for 
the odd-layered ones for both data sets. However, similar 
to what was previously found for the first data set, the 
coupling go ~ gi/N exhibits just the opposite size oscilla- 
tions. As Tc follows the trend ~ e~^/'^^'-°^\ one sees that 
the thickness-dependent oscillations of A^(0) are 7r-shifted 
with respect to the oscillations of 17, but the increase of 
A(0) has a more significant effect than the corresponding 
decreases of g for even number of MLs. It is interesting 
to note that the fourth and fifth experimental data-sets 
are characterized by almost the same coupling between 
electrons (the difference in g is less than 1%). So, for the 
experiments reported in Ref. [^, g is not sensitive to the 
change of Si(lll) surface crystal ordering from (7 x 7) 
to (-v/S X -s/S)- Here the question arises why the results 
of Ref. are so different with respect to those of Ref. [^. 
One possible reason may be the difference in the qual- 
ity of the silicon substrate which may have a significant 
impact on the interface phonons, as seen form the com- 
pletely different values of gi found from our analysis in 



(c) and (d). 

Note that the first and second data-sets in the present 
study were considered in our previous work.— However, 
in that paper we performed time consuming calculations 
invoking the full BdG formalism. As already explained 
above (see the discussion after Eq. ([9])), the full-BdG- 
based numerical procedure of determining the critical 
temperature often results in an overestimation of Tc that 
are larger than the errors of about a few percent induced 
by the Anderson approximation. This is reason why in 
the present paper we have chosen to use the Anderson- 
approximation-based analysis of the available data, in- 
cluding the two data-sets considered previously in Ref. 0. 

We also note that we did not discuss error bars for the 
experimental data in the present paper because they do 
not influence our results qualitatively: the fitting param- 
eters may be slightly different but the main conclusion 
about the crucial importance of the interplay of quantum- 
size oscillations of g and A^(0) will remain unaltered. 

IV. CONCLUSIONS 

We have analyzed five experimental data-sets for the 
critical temperature of ultrathin Pb nanofilms deposited 
on top of Si(lll). We have demonstrated that all these 
data sets can be reproduced within the BdG formalism. 
We addressed the important problem that the experimen- 
tally found amplitudes of the quantum-size oscillations of 
Tc vary significantly from one experiment to another. We 
showed that the reason for this difference is an interplay 
of the quantum-size variations in the single-electron den- 
sity of states with the thickness-dependent oscillations 
in the phonon-mediated coupling between the electrons. 
Such oscillations of the coupling depends on the substrate 
material, the quality of the interface, the presence or ab- 
sence of a protection cover and other fabrication details 
that change from one experiment to another. In addition, 
we demonstrated that oscillations in the coupling g and 
their interplay with the thickness-dependent variation of 
the density of states can even significantly change the 
overall trend of Tc with decreasing nanofilm thickness. 
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